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Abstract. Monte Carlo results for the moments (M k ) of the magnetization distribution of the nearest- 
neighbor Ising ferromagnet in a L d geometry, where L (4 < L < 22) is the linear dimension of a hypercubic 
lattice with periodic boundary conditions in d — 5 dimensions, are analyzed in the critical region and 
compared to a recent theory of Chen and Dohm (CD) [X.S. Chen and V. Dohm, Int. J. Mod. Phys. C 
(1998)]. We show that this finite-size scaling theory (formulated in terms of two scaling variables) can 
account for the longstanding discrepancies between Monte Carlo results and the so-called "lowest-mode" 
theory, which uses a single scaling variable tL d ^ 2 where t = T/T c — 1 is the temperature distance from 
the critical temperature, only to a very limited extent. While the CD theory gives a somewhat improved 
description of corrections to the "lowest-mode" results (to which the CD theory can easily be reduced 
in the limit t — > 0, L —* oo, tL d ^ 2 fixed) for the fourth-order cumulant, discrepancies are found for the 
susceptibility (L d {M 2 )). Reasons for these problems are briefly discussed. 

PACS. 05.70.Jk - 64.60.-i - 75.40.Mg 



1 Introduction 

Since about 15 years the five-dimensional Ising model is 
used as a "laboratory" [1-8] to test theoretical concepts 
about critical phenomena, in particular the concept of 
finite-size scaling [9-32], which has become an extremely 
valuable and indispensable tool for the study of phase 
transitions in condensed matter [22,33,34] and gauge the- 
ories of elementary particle physics [35, 36] . In this context, 
the d = 5 Ising model is of particular interest, since it ex- 
ceeds the upper critical dimension, d* = 4, and hence the 
Landau mean-field exponents exactly describe the critical 
behavior [37,38]. Also correction-to-scaling exponents [39] 
are known precisely [38] and fluctuations around the lead- 
ing mean-field description can be dealt with by simple 
perturbation theory; a renormalization-group treatment 
of fluctuation effects is not required here [37,38]. While 
for d > d* the hyperscaling relation dv = 7 + 2/3 (v, 
7, [3 being the standard critical exponents for correla- 
tion length £, susceptibility \ and order parameter M, 
respectively) does not hold and hence finite-size scaling in 
its standard version ("the linear dimension L scales with 
C [9-12,16,20,23]) does not hold either [12-14], a simple 
extension was proposed [1,2,17,18] which can be phrased 
as [2] "the linear dimension L scales asymptotically with 
a thermodynamic length £ T = (Af 2 /xb) _1/d oc i" 2 / d ". 
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Moreover it was suggested that ratios of moments of the 
order parameter distribution, such as Q = (M 2 ) 2 / (M 4 ) or 
(|M|) 2 /(Af 2 ) can easily be calculated from the so-called 
lowest-mode approximation [17], which was believed to be 
exact for the limit L — > 00, t — > 0, L/It fixed and should 
yield, apart from scale factors, universal finite-size scaling 
functions of L/£t- 

In view of these rather definite predictions [17], appar- 
ent discrepancies between the theoretical results and the 
Monte Carlo simulations [1-3] have been disturbing and 
it has been a matter of debate [27,4-8,29-31] whether the 
discrepancies reflect corrections to finite-size scaling. In 
Rcf. [8] it was shown that the Monte Carlo data for Q 
are indeed compatible with the predictions of Brczin and 
Zinn-Justin [17] if one takes into account two, theoreti- 
cally predicted, corrections to scaling. However, this still 
left the very slow convergence of Q as a function of L to- 
ward its predicted asymptotic value as a remarkable fea- 
ture (in Ref. [3] the data for Q for the considered range 
of system sizes appeared to have a common intersection 
point unequal to this value). More importantly, the cor- 
rectness of the treatment in Ref. [17] has recently been 
doubted [31] (see below). This controversy is cumbersome 
because the fact that it is apparently very difficult to dis- 
entangle the leading and subleading terms in finite-size 
scaling analyses even in a case where all involved criti- 
cal exponents are known precisely naturally leads to some 
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doubt on analyses where one wants to extract unknown 
critical exponents from finite-size scaling [22,23,26,34]. In 
addition, the problem also is of interest in the context 
of physical systems that are nearly described by Landau 
theory, such as systems with a long but finite range of 
interaction [25,28], polymer mixtures near their critical 
point of unmixing [40], etc. 

New light has been shed on this state of affairs by 
Chen and Dohm (CD) [29-31], who presented detailed ar- 
guments that for d> d* the standard treatment of the <j) 4 
field theory in continuous space [17,18,38] yields a mis- 
leading description of finite-size behavior, different from 
the finite-size behavior of a 4 model on a lattice, which 
one believes to be equivalent to an Ising model [37,38]. 
Chen and Dohm emphasized that therefore the justifica- 
tion given for the lowest-mode theory is invalid, and stated 
that the moment ratios mentioned above "do not have the 
universal properties predicted previously and that recent 
analyses of Monte Carlo results for the five-dimensional 
Ising model are not conclusive" [31]. 

In view of this criticism, a reanalysis of the available 
Monte Carlo results (including also some recent unpub- 
lished results [41] used in [42]) is clearly warranted. Such 
an analysis, where we compare the Monte Carlo data in 
detail with the result of the CD theory (which treats or- 
der parameter fluctuations perturbatively to one-loop or- 
der for the Ising case [31]) is given here. For the sake 
of a coherent presentation, we summarize the pertinent 
theoretical results in Sec. 2, while Sec. 3 gives a detailed 
comparison of the results for Q and for the susceptibility 
X = L d (M 2 } with the CD theory. Section 4 summarizes 
our conclusions. 



2 Theoretical background 

The singular part of the free-energy density //, of a system 
with linear size L in an external field h is written as (see, 
e.g., Ref. [14]) 



and 



f L = L- d f (tL Vt , hV> b , uL v ' ) , L —> oo 



(1) 



where t — T/T c — 1 and u is an irrelevant variable, in the 
renormalization-group sense (exponents y t > 0, j/h > 0, 
Hi < 0). Now for d > d* = 4 we have y t = 2, y h = (d+2)/2 
and yi = 4 — d, but u is a "dangerous irrelevant variable" 
(see, e.g, Refs. [43,44]), which means that the scaling func- 
tion f(xi,X2,Xs) is singular in the limit x 3 — > and can- 
not simply be replaced by f(x\, X2, 0). In terms of the bulk 
correlation length £b = Co^ ^ (above T c in zero field), the 
first argument of Eq. (1) can be interpreted as (L/^) 2 . 
Taking suitable derivatives of Eq. (1) with respect to the 
field we can thus write for the order parameter, the sus- 
ceptibility and the ratio Q (in zero field) 



(|M|) = L-^' 2 P M {t(L/£ ) 2 , (L/l ) A - d ) 
L d (M 2 ) 



d 2 f L 



dh 2 



= L 2 P x {t(L/to) 2 ,(L/£ ) 4 - d } 



Q = 



(M 2 ) 2 



P Q {t(L/to) 2 ,(L/£ ) 4 - d } 



(4) 



where Pm, P x an d Pq are the (universal [22]) finite-size 
scaling functions of the quantities (|M|), X and Q- The 
correlation-length amplitude £ of the bulk correlation 
length now appears as a scale factor for the variable X\ 
and for the variable X3 we have introduced the correlation- 
length amplitude £ of the bulk correlation length at T c 
in a small field [31] as a scale factor. In this way, the 
arguments of the scaling functions Pm, P x and Pq are 
dimensionless, as they should be. Note that u cx £q^ 4 - 

For large L the variable x 3 cx (L/£ ) 4 ~ d clearly be- 
comes very small, and hence it is an obvious question to 
ask how all these functions /l, Pm, P x , Pq behave in 
the limit x 3 — > 0. It was assumed in Ref. [1] that then 
the dangerous irrelevant variable X3 enters in the form of 
multiplicative singular powers of X3, e.g., 



f L (xi,X2,x 3 ) =xl 3 f L (xixl 1 ,x 2 x^ 2 ) . 



(5) 



This assumption was in the first place motivated by the 
fact that this is the mechanism that operates for the scal- 
ing in the bulk for d > 4 [44], and secondly by various 
phcnomcnological arguments. In particular, it was argued 
that standard thermodynamic fluctuation theory requires 
for T < T c and sufficiently large L that the distribution 
function Pj, (M) of the magnetization per spin for M near 
the spontaneous magnetization ±Mb is a sum of two Gaus- 
sians [1,2,12], 



Pl{M) = 



L d/2 



{exp [-(M-M b ) 2 L d /2 Xh \ 



+ exp[-(M + M b ) 2 L d /2xb]} . 



(6) 



Using M b - M h \t\? = Mt(-i) 1 / 2 and Xb = x b l^ 7 = 
X' b (— t) 1 the arguments of the exponential functions have 
the form 



(M/M b 



-1/2 



T-l 



(L/t T y 



(7) 



M b 2 ~x'ht 2 ■ Taking moments of this 



1 

2 

with £ d = M~ 2 Xb 
distribution one hence expects that the scaling functions 
in Eqs. (2)-(4) reduce to scaling functions of a single vari- 
able (L/£t)^ 2 oc tL d l 2 or, keeping the amplitudes £0 and 
£ explicitly present, tL d ' 2 ^ 2 £ (4 ~ d)/2 , i.e., 



|M|) = L-^Pm (tL d /% 2 £ 



-2M-d)/2 



'0 



) 



(2) and 



x = L d ' 2 P x (tL d /% 2 £t d)/2 ) 
Q = P Q (tL d /% 2 £t d)/2 ) ■ 



(8) 
(9) 

(10) 



Note that scale factors for the magnetization and the sus- 
(3) ccptibility have been absorbed in Pm (or Pm) and P x (or 
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P x ), respectively, while in ratios such as Q (and hence in 
Pq and Pq) such scale factors are divided out and fully 
universal functions remain. 

These arguments as they were presented in Refs. [1, 
2] did not tell anything about the explicit form of the 
scaling functions Pm(x), P x {x) an d Pq(x), however, and 
hence no prediction for the supposedly universal constant 
Pq(0) was made. In fact, making the (premature!) as- 
sumption that linear dimensions L = 3 to 7 lattice spac- 
ings are already large enough to obtain the limit x% oc 
(L/^o) 1 —> in d = 5 dimensions, it was argued that at 
T c there occurs a distribution of the scaled order param- 
eter Pl{4>) oc exp(Acj) 2 — 4 ), which implies a shift of T c 
as T c (L)/T c (oo) - 1 oc AL~ d l 2 , if T C (L) is defined as the 
temperature where Pl(<A) starts to develop a two-phase 
structure. However, the next step in the development, due 
to Brezin and Zinn- Justin [17], suggested that in the scal- 
ing limit Pl(4>) oc ex P(~ <t> A ) at T c , since the shift of T c as 
denned above should only exhibit a scaling with a higher 
power of L _1 , namely T C (L)/T c (oo) — 1 oc L 2 ~ d , because 
it results from corrections to the scaling description given 
in Eqs. (8)-(10). According to Ref. [17], the asymptotic 
behavior is simply given by the homogeneous order pa- 
rameter M in the finite system, 



P L (M) = exp 



-L d Qr-oM 2 + uM' 



(11) 



where r = a t changes sign at T c , u is the dangerous 
irrelevant variable mentioned above, both ao and u being 
nonuniversal constants. From Eq. (11) it is straightforward 
to derive that [17] 



P Q (0) = 87r 2 /T 4 (l/4) « 0.456947 . 



(12) 



However, the statements of CD imply that the continuum 
model considered in [17] actually leads to a behavior dif- 
fering from Eq. (11) and thus at this point also Eq. (12) 
seems without safe foundation. CD obtain, in the large-n 
limit of the n-vector model on the lattice rather than in 
the continuum, a result for the scaling function P x (x,y) 
[Eq. (3); x = t(L/£ ) 2 , y = {L/t Q ) i - d ] which is believed 
to be asymptotically exact, namely 



Px(x,v) = -r 6{x,y) + ^/[S(x,y)} 2 +4y 

JO L 



■i -1 



(13) 



where Jo is the interaction range of the </> 4 model on the d- 
dimensional hypercubic lattice (the lattice spacing being 
taken as unity here throughout), 



dL d 



and S(x, y) is given by 

S(x,y)=x-yh(J - 1 p- 1 ) , 
with the function I m (x), m= 1,2,..., being 
1 



(14) 



(15) 



I m (x) 



(27T) 



2 m 



/ 

Jo 



dyy m 1 exp(-xy/4:Tr 2 ) 



d/2 




(16) 



In terms of the Hamiltonian of the n-vector model with 
n-component vectors on the lattice, 

H = E [y^ + u °(tif] + E J f - ^) 2 . ( 17 ) 

the characteristic lengths ^o, m Eqs. (2)-(4) are given 

by 



e = ^d+5 c b ), c 4 - 4uon 1 



a 



J 2 l + 5 c b 



, (18) 



with r = r 0c + ao* and 

= u n J dk [8J(k)}- 2 , 5J(k) = J(0) - J(k) , (19) 

where J(k) = L~ d J^i j Jij ^exp[— ik ■ (rj — r^)]. For the 
n-vector model with n = 1, which is supposed to belong 
to the Ising universality class, comparable results are ob- 
tained only to one- loop order in perturbation theory [31]. 
Although the results are not exact, their scaling structure 
is analogous to Eqs. (13)-(16) and this structure is not 
expected to be changed by the higher-order terms of the 
loop expansion. Defining reduced moments 



CD find [31,45] 



J °° d<f><j) m exp [~\Y<j) 2 - <j) 4 ] 
/ °°# exp[-iY>2_ 04 ] 



Pq(x,v) 



1 6 2 (Y(x,y)) 
J o ^y + 36I 2 (f)y 2 
[02(Y(x,y))} 2 
9 4 (Y(x,y)) ' 



(20) 

(21) 
(22) 



with 



Y(x,y) =\x- 12y/i(f) - 1440 2 (a;y- 1 / 2 )/ 2 (f)y 3 / 2 / 



[y + 36I 2 (f)y 



2ll/2 



(23) 



where f = x + I202(xy~ 1 / 2 )y 1 / 2 . As should be clear from 
what has been said above, the results (21)-(23) should 
hold for sufficiently large L. 

Armed with these results we are now in a better posi- 
tion to reconsider the question already posed in Rcf. [1], 
namely to take the limit y — > 0. For this purpose we first 
consider the large-n limit, where we can rewrite Eq. (13) 
as 



i 



8(%y) 

2V? 



+ Jl + {6(x,y)/(2^yW 



(24) 

In the limit y — > we see from Eq. (15) that P x (x,y) 
depends, apart from the prefactor, on x and y through 
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the variable 

S(x,y) x 



2jy 2^y 2 



y/VllV^P- 1 ) 



= \tL*l%Ht d)/2 \{L/l^-^h{J^P^) . 

(25) 

Thus we see that there exists a limit t — * 0, L — > oo, 
tL d / 2 fixed, where the susceptibility scales exactly as pos- 
tulated in Eq. (9), since then the correction term of or- 
der (L/£o) {i - d)/2 in Eq. (25) clearly is negligible (J^P" 1 
tends toward zero in this limit, so the function I\ ap- 
proaches a finite constant). Contrary to statements made 
by CD themselves, viz. that the structure of Eqs. (8)- 
(10) is incorrect for the <fi 4 lattice model, we rather think 
that they have proven(!) the correctness of Eq. (11), in 
the limit specified above, at least for the large-n limit, 
and gratifyingly there is no contradiction at all between 
Eqs. (13)-(19) and the ideas of Rcf. [1] that led to Eqs. 
(8)-(10). Of course, the strong merit of the CD treatment 
is that it yields not only the scaling structure but also 
the explicit scaling function and a full description of the 
corrections due to the dangerous irrelevant variable u - 

We arrive at similar conclusions in the case n = 1, 
though one must recall that these results are only based 
on a one-loop order approximation. In the considered limit 
y — ► the quantity Y(x, y) in Eq. (23) reduces to 



Y(x,y)^—[l-18I 2 (f 



-12Ji(f)V» 

tL d '%H { t d)/2 [l-18/ 2 (f)(L/4) 4 - d ] 
-12h(f)(L/£ )^/ 2 , (26) 



which is, apart from the additional 0(L ) correction, 
fully analogous to Eq. (25). In the limit of interest (t — ► 0, 
L — > oo, tL d ' 2 fixed), f vanishes and the functions I\, I 2 
take finite values, so we see again that Eqs. (9) and (10) 
hold. Moreover, one concludes that at the critical temper- 
ature Y(0,y — > 0) — > and hence also Eq. (12) holds, as 
noted already by CD. It remains to be seen whether (12), 
which is less general than the scaling structure of Eqs. 
(8)-(10), holds to all orders in the loop expansion. 

In order to compare Eqs. (20)-(23) to numerical Monte 
Carlo data, it is clearly of interest to consider simple lim- 
iting cases of the susceptibility, where then the nonuniver- 
sal parameters £ an d £0 can be extracted. Since accurate 
Monte Carlo estimations of correlation lengths are much 
more difficult to perform than estimations of the suscep- 
tibility we note that in the large-n limit [Eqs. (13)-(19)] 
the bulk susceptibility is [29] 



Xb 



d 2 f h (t,h) 

dh 2 



Ik 

Jo 



1 + Sc b 

a t 



(27) 



where £b = £ot~ v - The same result is obtained from Eqs. 
(3), (13) using that, at fixed small t, d(x,y) w x in the 
limit L — > 00 and hence 



P x (x, y - 0) « (Jox)- 1 => X = L 2 /[J {L/i h f 



Xb 



In contrast, at the critical temperature the result is 

Jox(T = T c ) = L d ' 2 l { t d) ' 2 ■ (29) 

Thus, one can determine both parameters of interest, £0 
and £0, from the behavior of \ m two simple limits. The 
same procedure can also be carried out in the n = 1 case, 
considering the limit y — > at fixed small t, 



xt = 



e 



Jo v /y + 36/ 2 (f)y 2 



e 2 (Y(x,y)) y -^$_, 



Jo 



(30) 



while in the finite-size scaling limit (x = 0, y small) one 
obtains for d — 5 



L d i 2 r(|) 
JoVr r^) 



(31) 



(28) 



As we already noted, the interest of Eqs. (20)-(23) is not 
only the combination of the scaling structure, Eqs. (8)- 
(10), in the limit t — > 0, L — > 00, tL d l 2 = const, but these 
equations also incorporate the effect of the corrections to 
the lowest-mode approximation, which we would recover 
if in Eq. (23) we had Y(x, y) = xy' 1 ' 2 . 



3 Comparison of the Chen-Dohm predictions 
with simulation results for the d = 5 Ising 
model 

In Rcf. [8] only numerical data for the amplitude ratio Q 
have been considered, with 5 < L < 22. In order to es- 
timate the scaling parameter £o we now analyze the cor- 
responding data for the magnetic susceptibility. Thus, we 
apply a finite-size expansion similar to Eq. (3) in [8], 

X (T, L) = L d ' 2 (p + Pl tL^ + P2 i 2 L 2 y"t 

+ qi L*+q2L 2 *) , (32) 

where t = t + aL Vi ~ Vt and y* = y t — y\/2. So the term 
tL y t = tL d l 2 + aL^ d ^l 2 just corresponds to the scaling 
variable in (25) and (26). The additional term in Eq. (26) 
was already mentioned in [8] as the "cross term" tL Vt +Vi ; 
in contrast to the analysis of Q, it turns out to have a neg- 
ligibly small coefficient in the analysis of the susceptibil- 
ity. The leading power L d / 2 [Eq. (31)] has been confirmed 
numerically within less than one percent in Ref. [42]. In 
our analysis we have kept this power as well as the irrel- 
evant exponent y\ fixed. This yielded a critical coupling 
J/k^T c = 0.1139152 (4), in excellent agreement with the 
value found in Ref. [8] from keeping Q fixed at the zero- 
mode prediction (12), viz. J/k B T c = 0.1139150 (4) (num- 
bers in parentheses denote the uncertainty in the last dec- 
imal places). Furthermore, we found y t * = 2.53 (4), very 
close to d/2, and po = 1.86 (7). The quality of the fit in 
terms of the x 2 criterion was % 2 /DOF = 1.06. In order to 
improve the accuracy of our estimate for po, we have re- 
peated the analysis with y% fixed at d/2, finding J/k B T c = 
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Fig. 1. Log-log plot of l/(tx) versus tL 
dimensional nearest-neighbor Ising model 
Monte Carlo data as mentioned in the text, dashed lines rep- 
resent Eqs. (30) and (31), respectively, with the parameters 
from Eq. (33). The finite-size asymptote included here refers 
to the limit L — ► oo and has been estimated as described after 
Eq. (32). 



0.1139155 (2) and p = 1.91 (2) (x 2 /DOF = 1.05). All 
analyses were obtained with 5 < L < 22; upon omitting 
the smallest system sizes, a very similar estimate for po 
was obtained, with a minor increase in the uncertainty. 
For a more detailed analysis we refer to [41,42]. From our 
estimate for p and Eq. (31) we find, using J = 2J/k^,T, 
£o = 0.603 (13). For the sake of clarity, it is stressed that 
this estimate for £q thus pertains to the thermodynamic 
limit and is not a finite-size quantity. 

It is also possible to extract £o from the Monte Carlo 
data. However, here we use the series-expansion result 
from Ref. [46] for this purpose. Assuming the mean-field 
exponent 7 = 1, it was found that asymptotically \ = 
A/(l - v/v c ) with A = 1.311 (9) and v = tanh(J/fc B T). 
Rewriting this in terms of the reduced temperature t, we 
have X = A[cosh(J/fc B T 1 c)sinh(J/fcBT c )/(J/A : Br c )]t- 1 = 
1.322£- 1 and Eq. (30) shows that £ = 0.549 (2). 

Figure 1 now shows the log-log plot of (ix)" 1 versus 
the scaling variable tL d l 2 , using data for T > T c only 
(in view of the very accurate estimates of the critical cou- 
pling, the errors due to the inaccuracy of T c are not of ma- 
jor concern here). Available data for smaller system sizes 
have been omitted from this graph, because the rather 
strong deviations from scaling noted already in Ref. [2] 
would obscure its main purpose, namely to illustrate the 
use of the limits (30), (31) to extract £0 and £o- Of course, 
due to the corrections to scaling included in Eq. (32) the 
Monte Carlo data for L = 12 should not converge to the 
finite-size asymptote for L — > 00, but to a slightly shifted 
straight line. However, on the scale of Fig. 1 the finite-size 
asymptotes for L = 12 and L — > 00 arc indistinguishable. 
Because of their central interest we repeat our estimates 
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(33) 



Fig. 2. Plot of the moment ratio Q = (M 2 ) 2 /(M 4 ) for d = 5 as 
a function of the two variables x — tL 2 /£q and y = (lo/L) d ~ 4 , 
according to Eqs. (20), (22) and (23). 



The amplitude Qy/Io which normalizes the scaling vari- 
able tL d / 2 [cf. Eq. (26)] becomes 0.234 (4). 

In the following graphs, also Monte Carlo data from 
Ref. [2] (L = 4) and Ref. [3] (L = 8, 12) arc included and 
it was found that all the Monte Carlo data are, within their 
statistical errors, nicely compatible with each other. We 
have omitted the data of Ref. [3] for L > 13 here, since the 
rather irregular behavior found for the specific heat and 
the cumulant intersections for these system sizes indicates 
that these data suffer from statistical inaccuracies due to 
critical slowing down. Note that Ref. [3] used a singlc- 
spin-fiip Metropolis algorithm, whereas in Refs. [8,41,42] 
a single-cluster algorithm was applied. Available data from 
Refs. [4, 5] are restricted to temperatures very close to T = 
T c and hence arc unsuitable for our purposes. 

We now focus on the quantity Q, Eq. (4), since the 
scaling behavior of this quantity has been so controversial. 
Figure 2 gives a plot of the CD function (22), keeping both 
x and y as separate variables. One can see that for x large 
and negative Q = 1 as it must be and for x large and 
positive Q = 1/3, irrespective of y. This simply reflects 
the trivial properties of the low- and high-temperature 
phases, respectively. For \t(L/£ ) 2 \ < 1, however, a clear 
y dependence is seen. 

In Fig. 3 we compare the various Monte Carlo data to 
the CD function for Q as a function of tL 2 (i.e., the vari- 
able proportional to x). Note that in these plot there are 
no adjustable parameters whatsoever, so the agreement is 
at first sight very remarkable. At second sight, however, 
one does observe that there are slight but systematic de- 
viations between theory and simulation, which have con- 
sequences for the intersection of the amplitude ratios for 
different system sizes. Figure 4 demonstrates that both 
the CD function and the Monte Carlo data behave qual- 
itatively similar: for a range of sizes (4 < L < 12) there 
is almost a common intersection point, but it occurs at a 
negative value of tL 2 and consequently the corresponding 
ordinate value Qint is significantly larger than the pre- 
dicted asymptotic value (12). While this spurious value of 
the CD function, for the range of system sizes considered 



6 



Erik Luijten et al: Finite-size scaling above the upper critical dimension revisited 




(a) 



CD 

Refs. [8,41,42] 
Ref. [3] 




(b) 




(c) 

Fig. 3. Plot of Q versus tL 2 for (a) L = 4, (b) L = 8 and 
(c) L = 12. The full curves denote the predictions of Ref. [31]. 
Monte Carlo data generated at specific temperatures, taken 
from Refs. [2] and [8,41,42] are shown as open or full squares, 
respectively, while the histogram extrapolation data of Ref. [3] 
are shown as a broken curve. 
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Fig. 4. Magnified plot of Q versus tL 2 near tL 2 = 0, cf. Fig. 3, 
to demonstrate the occurrence of spurious cumulant intersec- 
tions. 
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Fig. 5. Plot of Q as a function of tll^ 2 , including both Monte 
Carlo data for L — 4, 8, 12 and the CD functions. Note that 
the zero-mode curve, resulting from setting Y = xy^ 1 ^ 2 in Eqs. 
(20), (22) and (23), practically coincides with the CD curve for 
L — 12 already, since £o [Eq. (33)] is so small. 



here, is about Q- mt w 0.48, it lies around Q int w 0.52 
for the Monte Carlo data; Rickwardt et al. [3] quoted 
Qint ~ 0.49 (1), including data up to L = 17. The lesson to 
be learned from this graph is threefold: (i) One must not 
pay too much attention to the value of such a cumulant 
intersection if one does not have a sufficiently large range 
of linear dimensions at one's disposal, (ii) The CD func- 
tion is a nice analytical example of a function that does 
produce a spurious "intersection" , as pointed out already 
in Ref. [31]: Although it looks so convincing on the graph, 
one knows that in the asymptotic limit the intersection 
occurs at t = and yields Q w 0.457 [Eq. (12)]. (iii) The 
CD function produces the same trend as the Monte Carlo 
data only qualitatively, but not quantitatively. 

What is the consequence of these results for the asymp- 
totic scaling, Eq. (10)? This question is addressed in Fig. 5, 
where the data from Fig. 3 are replotted as a function 
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Fig. 6. Plot of Q(T C ,L) versus L. Data points are the Monte 
Carlo results from Ref. [8], including statistical errors. The 
broken curve is the empirical fit, as described in Ref. [8], while 
the horizontal line denotes the asymptotic result of Eq. (12). 
The CD prediction is shown as a dotted line. 

of the scaling variable tL d l 2 : It is seen that the data for 
L = 4 deviate from scaling in a systematic way, while for 
L = 8, 12 the data scale already rather nicely, although 
they arc still a little bit set off in comparison to the theo- 
retical scaling curves. Note that on these large scales one 
cannot distinguish the CD curve for L = 12 from the 
lowest- mode result! The general trend appears that in the 
neighborhood of T — T c Eq. (22) yields a too small value 
for Q(L). 

In order to highlight the differences, we now amend 
the plot of Q at T c as a function of L, which was shown in 
Ref. [8] , by the prediction that would follow from CD (see 
Fig. 6): One indeed observes that the result of CD under- 
estimates the differences between Q(T C ,L) and Q(T c ,oo) 
distinctly — it basically yields a 1/ \[L correction, while the 
additional 1/L term resulting from CD is very small, un- 
like the rather pronounced 1/L correction that was found 
in the the empirical fit of Ref. [8]. Another, more ten- 
tative, way to quantify the differences is by adjusting £q 
such that the CD curve yields a reasonable description of 
the numerical data. It turns out that a value as high as 
£o ~ 3.2 (instead of the estimated value 0.603±0.013) is re- 
quired to find some agreement in the region 12 < L < 22. 
It clearly must be waited for a loop expansion to second 
order — which will yield additional 1 /L corrections of so far 
unknown magnitude — before one can draw final conclu- 
sions about the agreement between theory and simulation 
(or lack thereof). 

While in the comparison of the temperature depen- 
dence of Q as predicted by CD theory to the simulation 
results we have seen most of the disagreement for t > 0, 
Fig. 7 shows that much more drastic deviations between 
theory and simulation occur for \ m the regime t < 0. The 
fact that for t > there is perfect agreement for L = 12 
is no surprise, of course, since these data have been of 
central relevance for the fit in Fig. 1 that yielded £o and 
£ - It is clear that perhaps a better overall fit of the data 
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Fig. 7. Plot of x versus tL 2 for (a) L = 4, (b) L = 8 and 
(c) L = 12. The full curves denote the predictions of Ref. [31], 
Eq. (21), while the squares are Monte Carlo data of Refs. [8, 
41,42]. 
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Fig. 8. Plot of xL~ 5 ^ 2 versus tll^ 2 including Monte Carlo 
data for L — 4 and L = 12. Broken curves shown corresponding 
CD predictions, full curve is the zero-mode result. 



Fig. 9. Plot of x{Tc)L~ 5/2 versus L (where J/fc B T c = 
0.1139150). The dashed curve is a fit to Eq. (32) as described 
in the text. The dotted curve is the CD result (21). 



is reached if one would fit £q to describe the behavior of 
X for tL 2 large and negative, but then the behavior for 
t > would deteriorate. Let us note in passing that the 
term "susceptibility" is just used for convenience here: Be- 
low T c the real (reduced) susceptibility is of course given 
byi d ((M 2 }-(|M| 2 )}. 

Figure 8 shows then a plot of x_L~ 5 / 2 versus tL 5 / 2 , 
comparing Monte Carlo data for L = 4 and L = 12 [8] 
with corresponding predictions of the CD theory and the 
"zero-mode" curve. Again systematic deviations between 
CD theory and simulations are apparent: while the the- 
ory [31] converges to the zero-mode result from above, 
the Monte Carlo results fall clearly below the zero-mode 
result and nearly coincide with it for L = 12. This coin- 
cidence can be understood from a closer consideration of 
x(T c )L -5 / 2 versus L (Fig. 9): After a rapid increase from 
below to a value already close to the asymptotic value, 
this quantity flattens around L = 12 and then slowly ap- 
proaches (not necessarily in a monotonic way) its limiting 
value. In the whole region shown, the CD curve qualita- 
tively disagrees with the data — this disagreement clearly 
cannot be remedied by a different adjustment of the pa- 
rameters, because a monotonic decrease (close to a 1/y/L 
behavior) is an intrinsic feature of Eqs. (20)-(23) and also 
occurs in the large-n limit [Eq. (24)]. The deviation at 
L = 22 cannot be explained from a mis-adjustment of 
£o, since both curves approach the same limiting value 
for L — ► oo, where all finite-size corrections must vanish. 
Thus, if £q would have been chosen such that the CD curve 
coincides with the Monte Carlo result for L = 22, a mis- 
match would have occurred at L — > oo, which is clearly 
impossible. 

Of course, discrepancies between finite-size data for 
very small linear dimensions (such as L — 4 and L = 8) 
and the CD theory [Eqs. (21)-(23)], which only fully cap- 
tures the leading zero-mode result and the first correction 
terms (of order L -1 / 2 ) to it, would not be an argument 
against the usefulness of the theory However, Figs. 6 and 9 
clearly reveal that even for L = 22 one is still far from the 



regime where the CD theory satisfactorily describes the 
MC data. 



4 Concluding remarks 

In this paper Monte Carlo results for five-dimensional 
Ising lattices have been reanalyzed and compared to re- 
cent theoretical predictions obtained by Chen and Dohm, 
in an attempt to clarify a somewhat controversial discus- 
sion. Our results can be summarized as follows: 

(i) The CD theory reduces in the limit t — > 0, L — > oo, 
tL d ' 2 fixed, to the scaling structure proposed originally by 
Binder et al. [1] and explicitly illustrates the mechanism 
of multiplicative "renormalization" of variables by a dan- 
gerous irrelevant variable. In addition, it yields both the 
asymptotic scaling functions of various moments of the 
order parameter distribution as a function of the variable 
tL d / 2 and the leading corrections to it, which are of order 
_L(4- d )/ 2 . However, the comparison with the Monte Carlo 
data indicates that strong subleading corrections (of order 
1/L for d = 5) are present as well, which are not predicted 
by the CD theory, and one would need much larger L than 
accessible here (L max = 22) in order that these sublead- 
ing corrections are visible. So a quantitative agreement 
between theory and simulation is still far out of reach! 

(ii) The question must be asked to what extent the 4> 4 
model on a lattice for d > 4 yields the same behavior as 
the Ising model. Of course, one can take the parameters 
ro — > — oo, uq —* oo in Eq. (17) in such a proportion 
that the model precisely reduces to the Ising model (for a 
discussion see, e.g., Ref. [47]). At this point, we have not 
attempted to deal with this problem. 

(iii) The CD theory yields a nice illustrative example 
how one can be misled by apparent cumulant intersections 
which converge to the exact result extremely slowly as 
L — > oo. It is rather likely that this is the reason for the 
difficulties noted in the Monte Carlo studies in Refs. [1- 
3]. Both the CD theory and the simulations give clear 
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evidence that for a full understanding of the problem a 
variation of parameters over a broad range is desirable, 
including the behavior both above and below T c , as well 
as at T c . Corrections to the leading scaling behavior need 
careful consideration, which was already one of the central 
messages of Rcfs. [27,6,8]. 
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